CAUCHY PROBLEM FOR THE BOLTZMANN-BGK MODEL NEAR A 

GLOBAL MAXWELLIAN 

SEOK-BAE YUN 

Abstract. In this paper, we are interested in the Cauchy problem for the Boltzmann- 
BGK model for a general class of collision frequencies. We prove that the Boltzmann-BGK 
model linearized around a global Maxwellian admits a unique global smooth solution if 
the initial perturbation is sufficiently small in a high order energy norm. We also establish 
an asymptotic decay estimate and uniform L 2 -stability for nonlinear perturbations. 
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1. Introduction 

The dynamics of a monatomic, non-ionized gaseous system is known to be governed by 
the celebrated Boltzmann equation. But the complicated structure of the collision operator 
has long been a major obstacle in developing efficient numerical methods [4J. In an effort 
to find a simplified model of the Boltzmann equation, Bhatnagar, Gross and Krook [2j, and 
independently Walender [28J, introduced the Boltzman-BGK model: 

t d t F + v-V x F = u(M(F)-F) 1 

(1.1) 

V ' F(x,v,0) = F (x,v), 

where F(x,v,t) for (x, v , t) £ T 3 x R 3 x R + is the particle distribution function representing 
the number density of particles in phase space at position x, velocity v and time t. T 3 
denotes the 3-dimensional torus R 3 /Z 3 . Ai(F) is the local Maxwellian defined as 

. ,/„w \ p(x,t) ( \v — U(x, t)\ 2 \ 
M(F)(x,v,t) = , ; -exp - l - — ^ /' , 
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where p, U and T denote the macroscopic fields constructed from velocity moments of the 
distribution function: 



p{x,t) = I F(x,v,t)dv, 
J«. d 

p(x,t)U(x,t) = / F(x,v,t)vdv, 

jR d 

3p(x,t)T(x,t) = / F(x,v,t)\v - U(x,t)\ 2 dv. 

Throughout this paper, we assume that the collision frequency v takes the following form: 

v = v r] ^T) = p r > T w , 

where we have suppressed the constant to be unity for simplicity. A wide class of non-trivial 
collision frequencies is encompassed by this model. For example, Aoki et al. pQ studied the 
collision frequency defined as 

^i,o = 9- 

On the other hand, the following model was considered in [6l [20], 152] : 

where u was chosen to be the exponent of the viscosity law of the gas. The constant collision 
frequency [H El El [HI [2H [22l EH [25] corresponds to 

^0,0 = 1- 

The relaxation operator is designed to share important features with Boltzmann collision 
operator. For exmaple, the relaxation operator satisfies the following cancelation property: 



v I dv = 0, 

2 



(1.2) / u(M(F)-F) ( v 

which implies the conservation laws of mass, total momentum and total energy: 

d 



Fdxdv = 0, 
d f 

(1.3) — / Fvdxdv = 0, 



dt 
d_ 
dt 



T 3 xl 

F\v\ 2 dxdv = 0. 

T 3 xR 3 



We also have the following celebrated H-theorem: 

(1.4) ^ I Flog Fdxdv = I v(M(F) - F) log Fdxdv < 0. 

From the numerical point of view, the BGK model considerably simplifies the situation 
in that it is sufficient to update the macroscopic fields in each time step. But mathematical 
analysis is not necessarily easier, because the relaxation operator involves more nonlinearity 
compared to the bilinear collision operator of the Boltzmann equation. In [21], Perthame 
et al. established the global existence of weak solutions for the BGK model with constant 
collision frequency. Regularity and uniqueness was then considered in |17[ [22] under the 
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local existence frame work. The result in [T7] was employed in [25] to prove the convergence 
in a weight L 1 norm of a semi-Lagrangian scheme developed in [U [24} [23l 126] . which is, as 
far as the author knows, the first result on strong convergence of a fully discretized scheme 
for nonlinear collisional kinetic equations. In near-a-global-Maxwellian regime, the global 
existence in the whole space M 3 was established in [3] employing Ukai's spectral analysis 
[30] . Chan [5] studied the global existence in torus using the nonlinear energy method 
developed by Liu, Yang and Yu [18] . In [5], however, the decay rate is not known, which is 
expected to be exponentially fast. 

The purpose of the present paper is two-fold: first, we obtain the well-posedness of the 
Boltzmann-BGK model near a global Maxwellian for a wide class of non-trivial collision 
frequencies. Secondly, we establish the asymptotic decay estimate and uniform L 2 -stability 
[151 116j . The main theoretical tool is the nonlinear energy method developed by Guo |12[ 
[TBI E] to investigate the well-posedness of various important collisional kinetic equations 
such as the Boltzmann equation or the Vlasov-Maxwell (Poisson)-Boltzmann equation. 

Brief comments on possible extensions of our results are in order. Our assumptions on 
collision frequency do not cover the velocity dependent models proposed in [271 ETJ , which 
involves additional technical difficulties. Cauchy problems for relaxation models describing 
ionized plasma also can be considered by extending the arguments of this paper. We leave 
these topics for future research [29] . 

Before we proceed further, we set some notational conventions here. 
• (-, •) denotes the standard L 2 inner product in T d x R rf . 





£2 and || • 1 1 £2 denotes I? norms in and x M% respectively. 




• Multi-indices a, (3 are defined by 



a 



[a ,ai,a 2 ,cx 3 ], (3 = [ft, fa, fa] 



and 



• The energy norm 



is defined as follows. 



Ill/Will = E HW)lk,> 



\a\ + \P\<N 



where N > 4. 
• We define the high order energy norm for / as 




where the constant v c is defined in Proposition 12.21 
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• Throughout this paper, C a ^... will denote a generic constant depending on a, b, ■ ■ ■ , 
but not on x, v, and t. 

The paper is organized as follows. In section 2, we investigate the linearization procedure 
of the Boltzmann-BGK model. In section 3, the main theorem is stated. In section 4, we 
present several important technical lemmas. Section 5 is devoted to establishing the local 
in time existence and uniqueness of smooth solutions. In section 6, we study the coercive 
property of the linearized relaxation operator. Finally, in section 7, we combine these results 
to obtain the global in time existence of the classical solution. 



2. Linearized BGK model 



In this section, we consider the linearization of the Boltzmann-BGK model around the 
normalized global Maxwellian: 

I ,.,a 



m(v) 



_ H 
-.e 2 



We first establish a technical lemma which will be frequently used in the sequel. 
Lemma 2.1. Let G be a function given by 

p\U\ 2 + 3 P T _ 3p_ 
\/6' 



G(x,t) = 



\/6 



Then the Jacobian matrix of the change of variable (p, U, T) 

( 

d(p, P U,G) 



(p, pU, G) is given by 



d(p,U,T) 



V 



and 



1 














Ul 


p 











U2 





p 








U3 








p 





\U\ 2 +3T-3 


2pUi 


2pU 2 


3pU 3 


3p 




V6 


V6 


V6 


V6 


( 1 











\ 



( d(p, P U,T) 
\d(p, P U,G) 



A 



1 
p 




1 
p 








1 



B\ F>2 B3 C J 



where 



A = 
Bi = 
C = 



2|[/| 2 -(|[/| 2 + 3T-3) 



1 + 



U1+U2+U3 
p 



3p - (\U\ 2 + 3T - 3) 
2pUj-{\U\ 2 + 3T -3) 
p(3p-(\U\ 2 + 3T-3)y 



3p- {\U\ 2 + 3T -3)' 



Proof. It can be verified by a straightforward, but very tedious and lengthy calculation. We 
omit the proof. □ 
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Before we proceed to the next lemma, we define an operator which plays an important 
role in the theory of kinetic equations: 

Definition 2.1. The macroscopic projection is defined by 

5 
i=l 

where {e^} is an orthonormal basis for five- dimensional linear space spanned by {\/m, v\^fm, 
V2\/m, V3y/m, \v\ 2 ^/m}: 

ei = y/m, 
e 2 = viy/m, 

(2.5) { e 3 = v 2 y/m, 

e 4 = V3y/m, 

M 2 -3 /— 



Proposition 2.1. Let F = m + \[mj '. Then the local Maxwellian Ai(F) can be linearized 
around a global Maxwellian m as follows 

M(F) = m + pfV^+ £ (l\Dl Pe ^ Ge) M{e)}{i-e) 2 de)(f,e i ){f,e j ). 

l<ij<3 J ° 

Here M(0) denotes 

Pft \y-ue\ 2 

where pq,Uq,Tq are defined by the following relations: 
pg = Op +(1-0)1, 

^ 2 ^ PeU e = OpU, 

Pe\Ue\ 2 + ^PeT e 3 Q f P \U\ 2 + 3pT 3 i 
2 2 pe = 6 \ 2 2 P \- 

Proof. We define f(0) as follows 

/(») a M(e( P , P u, + 

= M(pe,peU e ,G e ) 
= M(0). 

We note that / represents the transition from the global Maxwellian m to the local Maxwellian 
M: 

(2.7) /(l) = —!- e—^~ and /(0) = —==e~ 4" . 

7(2^TF M ^ 7(2^ 

We then apply Taylor's theorem around = to see 

(2.8) f(l) = /(0) + /'(0) + - #) 2 d#. 

JO 



SEOK-BAE YUN 



(i) /'(0): We have from Lemma 12. II and the chain rule 
/'(0) = ±M{6(p,pU,G) + (1 - 0)(1,O,O)) 

= {p-l.pa^D^u^M^pU.G) - (1 - 0)(1,O,O)) 



(p- 


1, P U,G)- 




PeU 
9,Pe 


9, Gq) \ 

T e ) ) 


-1 








/ 1 








\ 


( 1 \ 






1 










V\ 


(p- 







1 



















1 




v 3 






V o o 





& 

u 3 


/ 


v h ;- 3 / 



M 



(p — l)m + (pU)vm + G- 



V6 



fy/mdvjm+ y J fy/mdvjvm + 
Pfyjm. 



V6 



mdv 



\ \v\ 



V6 



-m 



(ii) 



/ /"($)(! ~~ 9) 2 d6: We have from the chain rule 
■/ o 



/"(<?) 



d 2 M 



(9(p-l,pU,G)-(l-9)(l,0,ty) 



(p - 1, pU, G){Df pg ^ Ge) M(9)}(p - 1, P U, Gf 



We then substitute (i) and (ii) into (|2.8p to obtain the desired result. 
We now consider the linearization of the collision frequency. 



□ 



Proposition 2.2. The collision frequency can be linearized around the normalized global 
Maxwellian as follows. 



v = v r 



pi 



where 



Vc = (l) and Up = 6 ^ / D 

Here v{9) denotes 



(pe,PeUe,G e ) 



v{9){\ 



m = pI iy, 

where pg,Ug,Tg are defined as in the previous proposition. 

Proof. Since the proof is almost identical to the previous one, We omit it. 



□ 



Since the exact form of D? pgUe G ^Ai(9) is too complicated to be written down and 
manipulated explicitly, we introduce generic notations which considerably simplifies the 



THE BOLTZMANN-BGK MODEL 



argument. We first note from the chain rule 

/ 



D 



(p, P U,G) 



M{9) 



V 



(pe,PgUg,G e ) 



^ + B 

Pe 



J_ 

Pe 

v-Uei _ U ei 
PeTe ~~pf 
v-Ue2 _ Ue2 
PeTe p'j 
v-Ue3 _ Ug 3 
peTg ~~pf 
v-Tg + c \v-Ug\' 2 -3Tg 



Tg 



|2_ 
2fT 



J 



/ d(pe,p 9 U e ,G e ) \ - 1 
\ d(p e ,Ue,T e ) J 



V 



(pg,Ug,Tg) 



M(6) 



(i) 



\ 



X7, , I v ~ u ei UMX 

V (pg,Ug,Tg)[-^r- ~ -pj 

V7 , J v~Ue2 Ug 2 

V \p0,Vg,Tg)\-J g -Te~ ~ ~pT 

V7 , ( v-Ue3 Ug 3 

V (pe,Ug,Tg){^ g -Tg- ~ ~pT 



\ V ' { Pe ,u e ,T e )yy 9 



+ B 



V-Tg 
Tg 



+ C 



\v-Ug\ 2 -3Tg 
27f 



M(0), 



J 



where A, B, C are rational functions of macroscopic fields defined in Lemma f2.ll Therefore, 
we can deduce from Lemma 12. II that there exist polynomials Pff, Rfj such that 



D 



(pg,pglfg,Gg) 



Ug,Ug,Tg) 



1,3 



Rtf{pe,T 9 ,g ) 



where 



3Tgp e 



v 



Ge = 3 + 3 Pe - \U\j 
, x n ) and Rfj l (xi, . . . , x n ) satisfy the following structural assumptions {%m)'- 

- v-.m-v P(f is a polynomial such that Pij(0, 0, • • • , 0) = 0. 
• ^Hjv{2) Rfj 1 is a monomial. 
More precisely, we have for a multi-index m = (mi, 7712, • • • , 1113) 



{U M l) P l f(x 1 ,...,x n ) 
(U M 2) Rtfix!,... 



V n r mi T 

mi ni2 
Q"m.Xi Xn 



rn-2 
2 



1 , where ao = 0, 



From now on, we assume that Pfi* and Rff are defined generically, which means the exact 
form may change from line to line. These generic notations simplify the calculation dras- 
tically and cause no problems if we only keep in mind the structural assumption {Km) at 
each step. We now simplify the notation further by defining 



Q 



M 

i,3 



P % f{pe,Ue,T e ,v-Ue) 



R%{ Pe ,T ,g e ) 



and 



QM 



1 kA \ v ~ T e\ 2 
Qffe 







"(1 



2 d6 



to see from Proposition 12.11 

(2.9) M(F) = m + V^Pf + ^Y. Q % </. e *> (/, e j 
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By the exactly same argument, we write the collision frequency as follows: 
(2.10) v = i>c + J2 (%(/,*), 

i 

where 



Jo 

= ^(/, e ,,) f D iP(hPeU(hGe) v(e)(i-e)de. 



We again introduce generic polynomials P^(x\j . . . , x n ) and R^(xi, . . . , x n ) such that: 

D { P e,PeU e ,G e )v(. 9 ) 



P^ P0 ,Ug,T e ) 



R»(p e ,T e ,g e ) 

and assume they satisfy the same structural assumptions. 

{H v l) : (xi, . . . , x n ) = } y dmX^x™ 2 ■ ■ ■ x™ n , where a / 0, 

(2.11) 

(U u 2) : Ri{x±, ... ,x n ) = Cx™ 
We summarize the argument so far in the following proposition. 
Proposition 2.3. The relaxation operator is linearized around m as follows. 

: (M (F) - F) = {y c + u v ) { (Pf -/)+E«MX^i)| 



m 

ho 



We now substitute the standard perturbation F = m + y/mf into (jl.ip and apply Propo- 
sition [27TJ [2?2] and [2]3] to obtain the perturbed Boltzmann-BGK model: 

^L + vV x f = Lf + Tf, 

f(x,v,0) = fo(x,v), (x,v,t) eT 3 xR 3 xR + , 
where fo(x,v) = F ^J^ . The linearized relaxation operator L is given by 

(2-12) Lf = v c (Pf-f), 

and the nonlinear perturbation T(f) is defined as follows: 

r(/) = Up Lf + (u c + u p ) Qf 4 J (f,e l )(f,e J ), 

= Tiif, f) - r 2 (/, /)"+ f 3 (/, /) + r 4 (/, /, /), 



where 



ri(/,s) = VpPf = Y l (X{f,ei)( g ,e j )e j , 
hi 

r 2 (/,5) = "pf = , £ l <X(f,e j )g, 
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U(f,g,h) = ^Q^Qf 4 k (f,e l )(g,e,)(h,e k ). 

First we note that the conservation laws f)2 . 13j) now take the following form: 

Lemma 2.2. Suppose f is a smooth solution of \2.12\) . Then f satisfies the following 
conservation laws. 

/ f\/mdxdv = 0, 

Jt 3 xK 3 



T 3 xJ 



(2.13) / fv^/mdxdv = 0, 

f\v\ y/mdxdv = 0. 



T 3 xR 3 

We close this section by recalling the following important properties of the linearized 
Boltzmann-BGK model. 

Lemma 2.3. The macroscopic projection 

5 

1=1 

is a compact operator from L 2 into L 2 . 

Proof. It follows directly from the fact that the kernel of each integral operator lies in 
L 2 (T 3 x R 3 ) □ 

Lemma 2.4. L satisfies the following coercivity property. 

(Lf,f) = -u c \\(I-P)f\\ 2 Llv 
Proof. By (Pf, (I - P)f) = 0, we have 

(Pf,f} = (Pf,Pf) = \\Pf\\l lv , 

which yields 

(LfJ) = u c (Pf,f)-u c \\f\\l lv 
= -v c \\{I ~ P)f\\\ v - 

□ 

3. Main result 

We are finally in a position to state our main result. 

Theorem 3.1. Let N > 4 and Fq = m + -^/m/o ^ 0. Suppose that /o satisfies the con- 
servation laws \2.13\) . Then there exist positive constants C , M , 5* and <5* such that if 
E(0) < M, then there exists a unique global solution f(x,v,t) to \2.12\) such that 

(1) The high order energy norm is uniformly bounded: 

E(t) < CE(0). 

(2) The perturbation decays exponentially fast: 

\\\f(t)\\\(t)<e~ SH ^/Em■ 
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(3) if f denotes another solution corresponding to initial data f satisfying the same 
assumptions, then we have the following uniform L 2 -stability estimate: 

\\f(t)-f(t)\\Ll v <e- M ||/o-/o||^. 

Remark 3.1. Extension of these results to collision frequencies of the following form is 
straightforward. 

i=o j 

4. Preliminary estimates 

In this section, we present several estimates on macrosopic fields which are crucial to 
develop the argument further. 

Lemma 4.1. Let \a\ > 1. Suppose E(t) is sufficiently small. Then we have the following 
upper and lower bounds for macroscopic fields: 

(1) l-^E(t)<p(x,t)<l + ^E(t), 

(2) \U(x,t)\<3y/E(ij, 

(3) \<T[x,t)<\, 

(4) \d a P (x,t)\<^/m, 

(5) \d a U(x,t)\ <C H E(t), 

(6) \d a T(x,t)\ <C H E(t), 

for some positive constant C\ a \. 

Proof. (1) We have from Holder inequality 

p = 1 + J fV^dv < 1 + ||/|| 2 < 1 + y/E^j. 

Similarly, we have 

P > 1 - J f^Mdv > 1 - ||/|| 2 > 1 - yfW)- 

(2) Since J mvdv = 0, we have by Holder inequality 

J(m + y/mf)vdv j fv^fmdv 

< \\\fh ^3 JW) P 

< 3y/E(t). 

(3) The estimate of T can be treated similarly as follows: 

f(m + ^h-f)\v\ 2 dv-p\U\ 2 
3p 

3 + J f\v\ 2 ^mdv 



< 



3p 

3 + 240^^311/112 



< 

3p 
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l + 80\/^ v / £(7) 
3 

£ v 

where we used the smallness assumption on E(t) and 



4 e-^dv < 240\/2^3. 



The lower bound can be estimated analogously as follows: 

/(m + y/rnf )H 2 cfo - p\U\ 2 
3p 

3- J f\v\ 2 ^Edv- p\U\ 2 



> 



3p 



> 3 - 15y/w y/E^j ~ (1 ~ y/W))\3VW) 



> 



3(1 - yfmfi) 

l-(5^F + 3) v / ^(t) 



1-y/W) 
1 

> - . 

~ 2 

We now turn to the derivatives of the macroscopic fields. 

(4) follows directly by the same argument as in (1) noting that 

d a p = d a (J m + f^dv^j = J d a f^dv. 

(5) We observe from U = Lhlp^a th a t 

i^i<q<4 E h 2i }{ E \9p\ 2i }---{ E \d a p\ 2i }p~ 2H - 

l<i<|a| 1<7<M 1<7<M 

We now employ (1) and (4) to see 

\d"u\ < , , , 

" (1 - v^(t)) 2H 

where we used for i > 1 

E\t) < E(t). 

(6) Similarly, we observe that 

\d a T\ < C H { E \P\ 21 }{ E \dP? l }-{ E ^P^} 

l<i<\a\ 1<7<M l<7<|a| 

x { E E i^i 2l }---{ E \ru\»} P -* M . 

l<i<\a\ 1<7<M 1<7<M 

This gives by (1), (2), (4) and (5) 



\d a T\ < 



(1-7^)) 2 H' 

□ 
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The following lemma can be proved in an almost identical manner. We omit the proof. 
Lemma 4.2. Let \a\ > 1. Suppose E(t) is sufficiently small. Then we have 



(1) 1 - y/Etf) < Pe (x,t) <l + y/E(f) 

(2) \U e (x,t)\<3^E(fj, 

(3) ~<T e (x,t)<^, 



(4) \d a Pe (x,t)\<^W), 

(5) \d a U e (x,t)\ <C\ a \E(t), 

(6) \d a T e (x,t)\ <C\ a] E(t), 

for some positive constant C\ a \. 

Having established the preceding estimates for the macroscopic fields, we can now prove 
the following crucial proposition for the nonlinear perturbation T(f). 

Proposition 4.1. Suppose E(t) is sufficiently small such that estimates in Lemma \4-1\ and 
Lemma \4-£\ are valid. Then we have 

(1) | J d$T(f,f,f)rdv\ <C ]T \\d^f\\ L ,J\d^f\\ L2 J\r\\ L , 

l«i| + l"2|<|a| 

+ C £ H^/llLiJI^/IUgllrlU, 

|ai|+l<*2|<M, 

+ C £ CII^VlUl.JI^/llLgll^/llLgllrllig, 

|ai|+[ct2[+|a!3| 

(2, . _ _ 

x 



(ri >2)3 (/,5)/) + (rWs,/)/) <c S up||5|| Li j|/ni 2 , 

X 

(T 4 (f,g,h)f) + (T 4 (g,f,h)f) + (T 4 (g,h,f)f) < Csup \\g\\ L 2 sup \\h\\ L i ||/||£ 2 _, 



(3) T 1>2 ,3(f,g)r + T h2 ,3(g,f)r < Csup \r\ sup \\f\\ L 2 \\g\\ L 2 

L x,v X,V X 

T4(/, 9, h)r + r 4 (g, f, h)r + T 4 (g, h, f)r < C sup |r| sup ||/|| L 2 sup \\g\\ L 2 \\h\\ L 2 v 

X,V x,v X X 

Remark 4.1. Note that, unlike the case of the Boltzmann equation, we need to impose the 
smallness condition on the high order energy to prove the estimates. 

Proof. (1) To prove (1), we should consider Ti(f) (1 < % < 4) separately. But for simplicity 
we only present the proof for T^f). Other estimates can be obtained in an almost identical 
manner. 

The estimate of T^f): We first prove the following claim: 
Claim: There exists a positive constant e = e(a, j3) such that 



\«-Vq\ 

e 2T b 



\y-Ver 1 M 
< Ce ( 2+E > T e 4 



THE BOLTZMANN-BGK MODEL 



13 



(Proof of the claim): 

We apply the differential operator d to Qffe 



\v-Ue\ 



r ka \ v - u e\ 2 I H 2 
0[Qffe 2T o + 4 



| u -e/ 9 | 



2 i. ,2 



2 | ,2 

_+M_ 
4 + 



4 to see 



2T fl "f 4 



+ 



2 ,,ja 

4 



I^(j, ,U 6 ,T 6 ,v-U 9 ,v) -^i! + M! 
— i e s 



R{f(pe,Uo,T e ,g e ) 



< 



P^( P e,Ue,T e , 1,1) 
Rff(pe,Ug,T0,g e ) 

\y-ue\ 2 i H 2 



where we used the upper and lower bounds of Lemma 14,21 with 

g e = 3 + 3p e -|[/ e | 2 -3T e 



> 3 + 3(1 - v^)) - (3/B(t)) S 



2 

> 1 



3A/^(t) - 9E(t) 



and 



|« - U e \ r e ( 2+£ > T e < C{(2 + e)T e } a < oo, 



i«-f e l 2 



| v | r e < 2 + £ ) T e < CV(|t; — l/©r + |l/©r)e (2+£)T » < oo. 

Then the induction argument gives the desired result. We now employ the claim and use 
Holder inequality to see 

|a?r 3 (/,/)r|<fo 
E 



< 



|ao|+|ai|+|c*2| 



(4.14) 



< C 



E 



|ao| + |ai|+|<*2| 
=\a\ 



c 



E 



T{Qfje + " 4 ^(ff^f^d^f^rdv 

\v-Ug\ 2 \v\ 2 

e (2+ e) T e + 4 (ff*if iei )(ff**f iej )rdv 
\\d ai f\\ L? \\0 a2 f\\L ? , I < ~rclr. 



|ao|+|ai[+[a2| 
=\a\ 

We apply Holder inequality again to see 

f \v-Ug\ . \v\ 2 

(4.15) / e * hdv < C 



\y-u$\ , H 2 

g (2+ £ )T fl ^ 4 



L 2 



MIl2 <c|H 



l 2 - 
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where we used 

\y-Ua\ | M 2 
e 2T e 4 



exp 
< C [ exp 



1 



2 + e)T e 2J 4-(2 + e)T e 



£4 



1 



L(2 + e)T e 2 J 4-(2 + e)T f 



(2 + g)r g 

2 - (2 + e)T e 
dv 



dv 



( , ;( - 4 -(2+ £W ) 3 <co 



In the last line, we used Lemma 14.21 We now substitute (|4.15p into (|4.14p to obtain the 
desired result. 

(2) In (1), we set a=/3=0. Then we have from (I4.14p 



/ 



Lf. \\9\\ L?, 



T 3 (f,g)fdxdv < C 

< C J \\f\\ Ll \\g\\ Ll \\f\\ Li dx 

< C S up|| ff || L? ||/||| 2 . 



\ v ~ u e\ i M 2 
e 4 fd v )dx 



We now take L 2 norms with respect to spatial variables to obtain the desired result. 
(3) Let 4> £ L 2 . Then we have from the same argument used in (1) 

<r 3 (/,<?)r,0> < C J \\f\\ Ll \\g\\ Ll \\rct>\\ Ll dx 

< Csup|r| / ||/||l2||3|| L 2 UW^dx 

x,v J 

< Csup|r| 

x,v 

Therefore, the duality argument gives 



h\\g\\hdx nth* 



\T(f,g)r\\ < Csup\r\^ J\\f\\ 2 L2 J\gf L2v dx 
< Csup|r|sup||/|| L 2|| ff || L 2 



□ 



5. Local existence 

In this section, we establish the local in time existence of classical solutions under the 
assumption that the high order energy E{t) is sufficiently small. This local solution will be 
extended to the global solution in the last section by combining the coercivity estimate of 
L and a refined energy estimate. 

Theorem 5.1. Let Fq = m + ^/rnfo > 0. Suppose /o satisfies the conservation laws A2.13\) . 
Then there exist Mq > , T t > 0, such that if T* < and E(0) < there is a unique 
solution f(x,v,t) to the Boltzmann-BGK \2.12\) such that 

(1) The high order energy E{t) is continuous in [0, T*) and uniformly bounded: 

sup E(t) < M . 

0<t<T* 
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(2) The distribution function remains positive in [0, T*): 

F(x, v,t) = m + \Jmf{x, v, t) > 0. 

(3) The conservation laws \2.13\) hold for all [0, T*]. 

Proof. We consider the following iteration sequence: 

d t F n+1 + v • V^F^ 1 = v n {M{F n ) - F n+1 ), 
(5.16) . , 

F n+1 (x,v,0) = F (x,v), 

which is equivalent to 

{dt + v V x . + v c }r +1 = v c Pf n + rcP), 
f n+1 (x,v,0) = f (x,v). 
Now the theorem follows easily once we establish the following lemma. 



Lemma 5.1. There exist M Q > and T* > suc/i tfiaf if E(f ) < ^ then E(f n (t)) < M 
implies E{f n+l {t)) < M for t G [0, T*] . 

Proof. We take <9^ derivatives of (|5.17|) to obtain 

{d t + v ■ V x + v c }d$f n+1 = - J2i d P v ■ V *}d a f n+1 + VcdpPd a f n + d%T(f n ) . 
(5.18) ' v ' ^ 

Hi 

We take the inner product of (|5.18p with d^f n+1 and estimate each term separately. 

(1) L: l.h.s can be calculated directly as follows: 

(2) i?i: We see from the following observation 

d Vi v ■ V x d a f n+1 = 8 Xi d a f n+1 

that 

(Ri,o$f) < Eii^ u -v,}a a r +1 ii^ji^r +1 iii4. 

< cinr +i (t)ni 2 

< C£ n+ i(t). 

(3) R 2 : We first note that 

\\dpPd«f\\ Llv <Cp\\d«f\\ Llv . 

Therefore, we have from Holder inequality and Young's inequality 

(R 2 ,d$f}< < c\\&»n\h xv + mr +1 \\li v 

< c-(iiiriii 2 + iiir +i in 2 ) 

< C(E n (t) + E n+1 (t)). 
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We now turn to the estimate of the nonlinear term. 
(4) R3: We have from Lemma 14.11 



(R*,9$f) < C f \\d ai nLlJ^f n h^f n+1 \\ Li dx 

I u, ,7lR3 
l"l|+|o2|<|a| 

KI + M<M/ R 
l/9a|<|/?| 

+ c y / 3 ii^ l /iu i ji^ 2 rik ? ii^ 3 riu ? ii^r +1 iu s ^ 

|ai|+|a2|+|et3| 
<l«l 

< c £ (Bup||a*rii^ +^11^^11^) /ji^rii^ii^r 4 - 1 ^ 

|ai|+|«2|<|a| 

< c y {snp\\d^n L 2+ S np\\d^ni 2 )(\\d^ni 2 +\\dp n+i \\i* 



|ai| + |a2|<|a| 

3 



< C(Ei(t) + El(t) + yfEjt)E n+1 (t) + E n (t)E n+1 (t)), 

where we assumed a\ to be the smallest index without loss of generality and used the 
following Sobolev embedding 

(5.19) H 2 (T 3 ) C L°°(T 3 ). 

We substitute all these ingredients into (|5.18p to obtain 

^ii^r +1 (*)ii! i ,„+-cii^r+ 1 (*)iili„ 

< c[E n+l {t) + {E n (t)f + {E n (t)f + (E n (t))* E n+1 (t) + E n (t)E n+1 (t)Y 

We then sum over a and /3 and integrate in time to see 
E n +i(t) < E n+1 (0) 

+ C J [E n+l {t) + {E n {t)y + {E n (t)f + [E n {t)Y E n+1 {t) + E n {t)E n+1 {t)}ds 
<M^ + C \t* sup E n+1 {t)+T* sup E n + T,( sup E n )* +T*( sup E n f 

* L 0<t<T t 0<t<T* 0<i<T» 0<i<T» 

+ Tj sup E n Y sup E n+1 (t) + Tj sup E n ) sup E n+1 {t)\, 
V 0<t<T* y o<r» V 0<i<T* ' 0<T„ J 

which yields 



(l - CT» - CT*/Mo - CT*M ) sup £ n+1 (t) < (~ + CT* + CT» + CT,M W . 



0<t<T» 

This gives the desired result for sufficiently small Mq and T*. 



□ 



We now go back to the proof of the theorem and let n — > 00 to establish the local in 
time existence of a smooth solution. To prove uniqueness, we assume that g is another local 
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solution corresponding to the same initial data /q. We then have 

{dt + v ■ V + i/ c }(/ - 5 ) = P(f -g) + r 1)2)3 (/ - 5 , /) + r 1)2)3 ( 5) / - g) 
{ ' ] +U(.f-9jJ)+r 4 (gJ-g,f)+T 4 (g,gJ-g). 
We recall from Proposition 14.11 and (|5.19p 

(ri,2,3(/ -g,f) + ^1,2,3(9, f -g),f -9) 

+(r 4 (/ - g, /, /) + r 4 ( 5 , f-gj) + U(g, g,f-g),f- g) 

± E (H^llil,, + W^Wl^ + 11*711*8,, + Wg\\ Llv )\\f-g\\l lv . 

|a|<2 

We now multiply f — g to both sides of (|5.20p . integrate with respect to x, v, t and use the 
above estimate to see 

\\f(t)-9(t)\\h + f !!/(*)- 



sup (jRtf) + Jl^) + E f (t)+E g (t) + l) f \\f(s)-g(s)\\ 2 Li ds. 



: f * ~ 

|a|<2 1 

Therefore, for sufficiently small Ef(0) and E g (0), the uniqueness follows from Grownwall's 
theorem. We now turn to the continuity of E(t). Let / be the smooth local solution 
constructed above: 

d t f + vVf + u c f = Pf + T(f). 
We multiply / and integrate over x, v and then over [s,t] to see. 

\E(t)-E(s)\<c[l + VW) + E(t)] f E WfUlJr^Q. 

The positivity of m + y/rnf can be verified iteratively from the positivity of -Fo using 
(|5.16p . Finally, since the local solution is smooth, the conservation laws can be obtained 
straightforwardly. □ 

6. COERCIVITY OF L 

The coercivity estimate in Lemma 12.41 involving only microscopic components, is not 
strong enough to play as a good term in the energy method. In this section, we show that 
the full coercivity L can be recovered as long as the energy E(t) remains sufficiently small. 
We first set for simplicity 



Recall that / can be divided into its hydrodynamic part Pf and microscopic part (I — P)f: 

(6.21) f = Pf + (I-P)f, 

where 

Pf = ay/rn + b ■ vyfm + c\v\ 2 \/m. 
We observe that there exists constants C such that 

( 6 - 22 ) ^W-P)f\\ L i v < W-P)f\\ L %, v < c\\{i-p)f\\ Llv . 
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Now, we substitute (|6.2ip into the BGK model (|2.12p to obtain 
(6.23) {d t + v V}Pf =£{(1- P)f} + h(f), 

where 

£{(I-P)f} = {-d t -v-V x + L}{I-P}f, 

Hf) = r(/). 

The l.h.s of ()6.23j) is calculated as follows 

{^cM 2 + (d t c + d l bi)v 2 + ^2(8% + djb i)^ v J + + dla ) v * + dta]V^, 

i j>i 

where d l = d Xi . We then expand the l.h.s of (|6.23p with respect to the basis: 

\/m, Vi\pm, ViVjy/m, v 2 \frn, Vi\v\ 2 ^Jm (1 < i,j < 3). 
Equating both sides of the above identity, we obtain the following result. 

Lemma 6.1. a, b, c satisfy the following relations. 

(1) Vc = £ c + h c , 

(2) d t c + d i bi=£ i + hi, 

(3) & i b j + dib i =£ ij + h ij , 

(4) d t b i + d i a = £ bi + h bi , 

(5) d t a = £ a + K, 

where i c , ti, £ij, lu, 4 o-re coefficients of the expansion of £ with respect to the preceding 
basis. Similarly, h c , hi, hij, h^, h a denotes the corresponding coefficients of the expansion 
ofh. 

For brevity, we define i and h as 

I = 4 + J> + J>i + ^4i + 4, 

i i,j i 

h = h c + ^ hj + hjj + y^ j h bi + h a . 

i i,j i 

Lemma 6.2. Let \a\ < N — 1. Then we have 

|a|<7V-l [a[<JV-l 

W \\d?bi\\ L} < £ \\d a i\\ Llv + £ \\8 a h\\ Llv . 

|a|<iV-l |a|<AT-l 

Proof. (1) Following [12j [T3l [H] , we observe that 

-Abi - didA 

= -Yd j {d j bi)-2d i d i bi 

= -Y. 8 j (-8% - tjj - hjj) - 2d\-8 t c + £j + hj) 

( by Lemma 16.11 (3) and (2)) 
= ^(d'd'bj + d%c) + ^(cVlv - d j h ij ) - 2d\£i + hi) 
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= ^2(di(-d t c -it- hi) + &dtc) + Y,( dj£ *J ~ Phi) ~ 2&{jk + hi 

( by Lemma 16.11 (2)) 
= ^2(di£i + d%) + ^(^'^ - - 2d l (£i + hi). 

Then the result follows from the standard elliptic estimate. 
(2) By Poincare inequality and Lemma 16.21 we have 

< \\V x d^b\\ L 2 

< C(||^i||^ + ||^^|| i2 J. 



Lemma 6.3. For |o| < N — 1, we have 

(1) l|c|Ui<C( \\I\\ Li + fhh* ). 

(2) \\d t d a c\\ Ll <C{\\d a £\\ Ll + \\d a h\\ Ll ). 

(3) \\V x d a c\\ Ll <C{ 113^11^ + 11^11^ ). 

Proof. (1) By Poincare inequality and Lemma l6.1f l). we have 

\\4 L 2 x < ||VC|| L 2 

< C\\£ c + h c \\ L 2 

< C[\\£\\ L 2 + fh\\ Ll ), 
where we used the conservation of energy: 



c(x)dx = 0. 

(2) By (2) in Lemma 16. 11 we have 

\\dtc\\ L 2 = || - V-b + £ + h\\ L 2 

< c(\\I\\ Li + fh\\ Ll ). 

(3) follows directly from (1). 

Lemma 6.4. Let \a\ < N — 1. For (2), we assume further that a is purely spatial: 
[0,01,02,03] +- 0. Then we have 

(1) \\d t d a a\\ Ll <C{ \\d a £\\ Ll + \\d a h\\ Ll ), 

(2) \\Vd%a\\ Ll < E H^IUS' 

|o|<AT-l |a|<AT-l 

(3) ||a|| Z 2 < \\£\\ Ll + \\h\\ L% . 

Proof. (1) This follows directly from Lemma 16. II (5). 
(2) By Lemma 16. II (4). we have 

Ad a a = V-Vd a a 

= V- (-d t d a b + d a £ b + d a h b ) 
= -Vd t d a b + Vd a £ b + Vd a h b . 

We then multiply Va to both sides and use integrate by parts to see 

||va>|| L? < Wd^bUz + H5-4II + \\d a h\\ Li 
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< lia^eP- 1 ^ + \\d a £ b \\ + \\d a h b \\ Li 

< E E ii^ii^.- 

\a\<N-l \a\<N-l 

In the last line, we employed Lemma 16.21 

(3) follows from Poincare inequality combined with (2) and the conservation of mass: 

llallra < II Vail ra 

Lemma 6.5. For \a\ < N — 1, we have 

(1) E \\vt\\n<c E WV-P^fhiy 

M<iV-l h\<N 

(2) E \\d a h\\ Ll <C{^W + M ) E \\d a fhl v . 

\a\<N \a\<N 

Proof. (1) Note that there exists constants \ n such that d a £ takes the following form 

13 



□ 



f]X n [ d a £{(I-P)f}-e n (v)dv, 



where e n denotes the orthogonal basis for the 13 dimensional space spanned by 

{y/m,Viy/m,vfy/m,ViVjy/m, \v\ 2 Viy/m | 1 < i,j < 3}. 
We then observe that 

|2 



Ll 



\\d a i({I-P}f)-e n (v)dv\ lL2x 

= | J ( - {d t + v • V + L}(I - P)d a f) ■ e n (v)dv 

< J \e n (v)\dvx 

J \e n (v)\{\(I - P)d°d a f\ 2 + \v\ 2 \(I - P)V x d a f\ 2 + \(L(I - P)d a f\ 2 ]dxdv 

< C{\\(I - P)d°d a f\\ Llv + ||(/ - P)Vaa/||^ + ||(/ - P)d a f\\ Llv } 2 . 

This completes the proof of (1). 
(2) As in (1), terms in d a h can be presented as 

13 „ 

EA n / d a T(fJ)-e n (v)dv 
for some constants A n We now apply Lemma 14.11 (3) to get 
d a T(f) ■ e n (v)dv 



E / Ti^(d ai f,d a *f)-e n (v)dv 

l l^l l 



I"i|+|a2|<|a| 

+ E 

|ai|+|aa|+|Q!3| 
<l«l 

< C{VM + M) \\d a fhl v - 

\a\<N 



Ll 



T A (d a 'f,d a2 f,d a2 f)-e n (v)dv 
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This completes the proof. □ 

We can now prove the main theorem of this section. 

Theorem 6.1. Let f be a classical solution of Then there exists M and 5 = 5(M) 

such that if 

Y \\d a m\\h < M ' 

\a\<N 

then There exists 5 > such that 

Y (Ld a f,d a f)<-5 ]T \\d a f\\l lv . 

\a\<N \a\<N 

Proof. By Lemma 16.21 - 16.4] we have 

Y \\d«a\\ Ll + H^IU, + \\d«c\\ Ll < Y, H^lk, + E W da hhly 

\a\<N |a|<AT-l \a\<N 

We then apply Lemma 16.51 to see 

E H^/ik, ^ c E + H^ik + W4» m 

\a\<N \a\<N 

< c Y \\(i-P)d a f\\ Llv + c^W Y \\d a f\\ Llv . 

\a\<N H<N 

Hence we have from Lemma 12.41 and the equivalence estimate (|6.22|) 
J2(Ld a f,d a f) = -u c \\(I-P)d a f\\l lv 

\a\<N \a\<N 

< -v cC y w-p)d a f\\i lv 

\a\<N 

< -v cCl { y wpft Llv -c 2 VM y ir/iiiu 

\a\<N \a\<N 

< - min{ "f cCll { E W9 a Pf\\h +\\d a d-P)f\\h } 

\a\<N 

+\u c c 2 Vm Y ll^/llia,, 

|a|<iV 

< _^{min{l,C7i}-C72>/M} E H^/HU,- 

|a|<7V 

We then choose M sufficiently small such that 

min{l,Ci} > C2VM 

to obtain the desired result. □ 



22 



SEOK-BAE YUN 



7. Proof of the main theorem 

In this section, we derive a refined energy estimate for the Boltzmann-BGK model and 
establish the main result. Let / be the unique smooth solution constructed in Theorem l5.ll 
First we take d a on both sides of (|2.12p to have 

[d t + v-v + L]d a f = d a r(f). 

We multiply d a f, integrate over T d x R rf and apply Lemma 14.11 (1) and Theorem 16.11 to 
obtain 



For j3 ^ 0, we take &% to obtain 

(7.24) [dt + V'V + v c ]d$f = -Y,d%+ B e if n+1 + dpPd a f + d%T(f), 

i 

where e x = (l,0,0),e 2 = (0,l,0),e 3 = (0,0,1) and ei = (0,l,0,0),e 2 = (0,0,l,0),e 3 
(0,0,0, 1). and we used the following relation: 



l<i<3 

i — r r\/v o X 

V 



i 

Multiplying daf to both sides of (|7.24p and integrating over T d x R d , we obtain by an 
almost identical manner as in the local existence case 

i 

+c E / II^/IIx^II^/IIliJI^/H^^ 

\ai\ + \a 2 \<\a\ 

+ c E j Rd II^/IIliJI^/IIliJI^/H^^x 

|ai|+|o2|<|a| 
Ifti|<l0l 

+ c E I d i ii^viuiji^/ii^ji^/iki.ji^/ii^/x. 

, , , , , , JK. d JM. d 



|ai| + |a2| + |«3| ' 



Therefore, we have for sufficiently small e 
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+C E / Jd ai fhlJd a2 f\\LlMfh lv dx 

|ai|+|a2|<|Qi| 

+ C E j Rd \\d ai f\\LiJd^f\\ Ll Jdp\\ Li Jx 

|ai|+|«2|<|a| 
\M<\P\ 

+c E / / ll^/IUiJia-zikiJI^/lkiJI^/ll^dx 

The estimates for / can be treated almost identically as in the proof of Theorem 15.11 to 
obtain 



I np <C{JW) + E{t)}\\\f\\\\ 
where we used the Sobolev embedding H s » L°°. This yields 

^ ^ii^/iik. + Tii^nk, * ^Eii^-ei/iili,,+^ii^/iili,, 

+ C{ x /^) + £(t)}|||/||| 2 . 
From the above inequality, we observe that the bad terms of J2\/3\= m +i E< 8i * na ^ ^ s 

E {^En^!/ii!^+ c £ \m\\i lv y 

|/3|=m+l « 

can be absorbed in the good terms of C m Y^\p\=m ^8 + C m Yla E a if Cm is sufficiently large. 
Therefore, by an induction argument, we can find constants C m and 5 m such that 

|a|+|0|<JV, 
|/3|<m 

We now suppose E < 1 without loss of generality and set m = N to obtain 

E {3^iwi£i„ + m^/ii£ 5 J^c^iv^Diii/iii 2 . 

\a\ + \P\<N 



Notice that we used E(t) < yjE{t) and redefined 2C N by C N . We now define y(t) as 



2/(0 - E ^N-^WiJWl, 



\a\+\P\<N 

We choose a constant Ci such that 

1 f w , 57V 



y{t) + irl \\\m\\\ 2 ^}< E (t)<ci{y(t) + s -fJ o \\\f( S )\\\ 2 ds} 



Ci 

We define 

M 
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and choose the initial data sufficiently small in the sense that 

E(0) <M< M . 

Let T > be given as 

T = sup ji : E(t) < 2Cf Af} > 0, 

which gives 

E(t) < 2CfM < M Q . 

We then have for < t < T 



y'(t) + 6 N \\\f\\\ 2 (t)<C N VW)\\\f\\\ 2 (t) 
(7.25) <CWCiV^M|||/||| 2 (t) 

< 5 -f\\\f\\\\t). 

Therefore, integration above over < t < T yields 

E(t) < Cl {y(t) + 6 -f J*\\\f(s)\\\ 2 ds} 

< C iy (Q) 

< CfE(0) 

< C\M 

< 2CfM. 

This is a contradiction considering the continuity of E and the definition of T. Hence we 
have T = oo. By (|7.25p and y(t) < C|||/(i)||| for some constant C, we have 

y'(t) + 6 -fy(t)<y'(t) + 6 -^\\\f\\\ 2 (t)<0, 
which gives the exponential decay of the perturbation. 

We are now left with the L 2 -stability estimate. Let / be another solution corresponding to 
initial data /o- We subtract the equation for / from the equation for / to see 

{d t + v-v}(f-f) = L(/-/) + r 1)2)3 (/-/,/) + r 1)2i3 (/,/-/) 

+ r 4 (/ - /, /, /) + r 4 (/ 5 / - /, /) + r 4 (/, /, / - /"). 

We then multiply f — f and integrate over x, v and t to have 

^4n/(*) - mwh = w - /),/- /-) + (r 1>2 , 3 (f - /,/) +r lj2)3 (/,/- /),/- /) 

2 at x < v 

+r 4 (/ - /, /, /) + r 4 (/, / - /, /) + r 4 (/, /", / - /), /-/). 

We now apply the coercivity estimate in Theorem 16.11 

(L(f-f),f-f)<-5\\f-f\\l lv 

and the following estimates from Proposition 14.11 and (15. 19f) : 

<ri, 2 ,3(/ - /, /) + r 1)2|3 (/, / - /"), / - /) 

+(r 4 (/ - /, /, /) + r 4 (/~, / - f, f) + r 4 (/, /", / - /), / - /> 

< c (\\d a f\\l lv + \\d a f-\\li v + + \\d a fhi J 11/ - f\\l lv 

l«l<2 

<c{y%(t) + ^E}(t)}\\f-f\\ 2 Llv 
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to get 

\j t \\m - mwi lv + (s- {y/E-,( t ) + y/Efi)}) \\f(t) - mwi lv < o, 

which gives the desired result for sufficiently small E g (0) and Ef(0). 
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